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We study the behaviour of inflationary density perturbations in the vicinity of horizon crossing, 
using numerical evolution of the relevant mode equations. We explore two specific scenarios. In one, 
inflation is temporarily ended because a portion of the potential is too steep to support inflation. 
We find that perturbations on super-horizon scales can be modified, usually leading to a large 
amplification, because of entropy perturbations in the scalar field. This leads to a broad feature 
in the power spectrum, and the slow-roll and Stewart-Lyth approximations, which assume the 
perturbations reach an asymptotic regime well outside the horizon, can fail by many orders of 
magnitude in this regime. In the second scenario we consider perturbations generated right at the 
end of inflation, which re-enter shortly after inflation ends — such perturbations can be relevant for 
primordial black hole formation. 
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I. INTRODUCTION 

The inflationary cosmology remains the leading can- 
didate theory for the origin of structure (see Refs. [0J|] 
for reviews). Considerable attention has been focussed 
on producing highly accurate calculations of the pertur- 
bations produced during inflation j3j, which ultimately 
may be measured at the percent level via the microwave 
anisotropies they induce Special interest has been 
given to the simplest case, where there is only a single 
scalar field degree of freedom during inflation, and we 
will focus exclusively on that case in this paper. 

In considering how a perturbation might evolve, a cru- 
cial quantity is the comparison of the inverse wavenum- 
ber with the Hubble length, given by the ratio aH/k. By 
definition inflation corresponds to a decreasing comoving 
Hubble length, d(aH)/dt > 0; consequently when modes 
of a given wavenumber are considered, they begin their 
evolution well inside the horizon^ and cross outside dur- 
ing inflation. At early times, flat space-time quantum 
field theory can be used to fix the initial normalization 
of the perturbations. Around the epoch of horizon cross- 
ing, the perturbation becomes frozen in, corresponding to 
a constant curvature perturbation once the mode is well 
outside the horizon. This leads to standard expressions 
for the spectrum of density perturbations such as 
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where Viz is the power spectrum of the curvature pertur- 
bation 1Z, using the notation of Refs. |IJ^]. This expres- 
sion gives the amplitude of the perturbations in terms 



of the values of the Hubble parameter H and scalar 
field velocity <fi at the time the mode crossed the hori- 
zon, i.e. when k = aH. A similar expression holds for 
the gravitational wave amplitude, and more sophisticated 
higher-order versions have also been derived. 

Eq. ([I]) gives the impression that the amplitude of per- 
turbations is being quoted at the instant of horizon cross- 
ing. However, it is important to realize that that is not in 
fact the case. The value quoted is the perturbation am- 
plitude attained in the asymptotic limit k j aH — > when 
the perturbation is well outside the horizon; it just hap- 
pens to be written in terms of the values the background 
parameters had at the instant of horizon crossing. In fact 
the value of the curvature perturbation at the instant of 
horizon crossing typically differs by a significant factor 
from the asymptotic value. 

Bearing that in mind, in this paper we consider two 
circumstances where the standard formula for the per- 
turbation amplitude may not be valid, due to a failure 
to reach the asymptotic limit. In each case, this is be- 
cause inflation ends before the true asymptotic regime 
has been reached. In our first scenario, we consider a 
temporary interruption to inflation, where the field driv- 
ing inflation has a region where the potential is too steep 
to sustain inflation. We will see that this can lead to 
significant modifications to the standard results, and in- 
deed that even modes significantly outside the horizon 
can receive a large change in amplitude due to entropy 
perturbations in the scalar field. In the second scenario, 
we consider perturbations produced at the end of infla- 
tion, where their amplitude on re-entry has relevance to 
the production of primordial black holes. 



II. FORMALISM 



*We use 'horizon' as shorthand for the more precise 'Hubble 
length'; it does not refer to the particle horizon. 



The scalar perturbations are best followed using the 
variable u = a Sip and the equation satisfied by its 
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Fourier modes uy_ is 
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where primes denote differentiation with respect to con- 
formal time, z = acj)/H, and 
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where \^ denotes differentiation with respect to cf>. Mode 
equation (||) has two asymptotic regimes characterized by 
the relative sizes of fc 2 and z"/z. Only in the slow-roll 
limit, where e,|77|,|£ 2 | <§; 1, will this necessarily be the 
same as comparing fc and aH . Bearing this in mind, well 
within the horizon in the limit of aH/k —>■ each mode 
behaves like a free field 
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while in the limit k 2 <C z" / z we have a growing mode 
solution 
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which means that the curvature perturbation, |7?-k| = 
|uk/z|, remains constant on superhorizon scales. The 
quantity z is sometimes described as the pump field for 
scalar perturbations. 

To calculate the constant of proportionality of Eq. (||) 
it is generally assumed that e and r\ are slowly varying 
at around horizon crossing, which will be a valid approx- 
imation as long as these parameters are small, since 
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The expression for conformal time r then takes on a sim- 
ple form 
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and Eq. (g) reduces to a Bessel equation leading to the 
Stewart-Lyth result for the power spectrum B 
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where C ~ —0.73 is a numerical constant. 

In fact, the constancy of IZ^ depends on e and r/ doing 
nothing dramatic even after horizon crossing, which can 
be seen if we rewrite Eq. (0) in terms of TZ^ itself 
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where 



-K' k + k 2 K k = . 



= aH[l + e- n\ 



(12) 



(13) 



Initially TZ k , like u kl will be oscillating, although we are 
only interested in the envelope of this oscillation, |72.k|i 
because only the envelope contributes to the value of 
the real space curvature perturbation, 1Z. Mode equa- 
tion (|l^) is of the form of a damped harmonic oscilla- 
tor. In the slow-roll limit, at around horizon crossing the 
system becomes dominated by the exponentially grow- 
ing friction term proportional to TZ k , and the solution to 
Eq. (Il2| ) soon becomes well approximated by the form 



TZk(T) = const 
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where TV is the number of e-folds after horizon crossing. 
The rapid freezing in of the curvature perturbation is ap- 
parent from Eq. (|l4|), which measures the rate of change 
of 1Z per Hubble time. We will examine the properties of 
the late-time solution in more detail later in the paper. 

In this paper we are interested in a more dramatic sit- 
uation, arising through a failure of slow-roll. If at any 
time after horizon crossing the friction term in Eq. (H3) 
changes sign and becomes a negative driving term, then 
we can expect dramatic effects on modes which have re- 
cently left the horizon. This change of sign will occur 
whenever z reaches a local maximum, or equivalently 
whenever 
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rj = 0. 
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As e is always positive, 77 must be at least one for this 
to happen, which implies that a turn around in z can 
occur only during a transition to fast-roll inflation or to 
a non-inflationary period. 

Before progressing to specific applications of the above, 
we note that no such interesting effects can occur for 
gravitational waves. Their mode equation is given by 
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In analogy with the above, Eq. jl^ ) can be written as 
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FIG. 1. The comoving Hubble wavenumber, aH, increases 
with the number of e-folds of expansion, N, during inflation. 
We set N = when inflation first ends, and inflation is sus- 
pended for around 1 e-fold. The quantity 1 + e — 77 remains 
negative for around 5 e-folds. We took B = 0.55. 



where 
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The pump field a increases for all time, and so the con- 
stancy of the gravitational wave amplitude after horizon 
exit is assured, until horizon re-entry. 



III. TEMPORARY INTERRUPTION OF 
INFLATION 

To illustrate the points made in the previous section, 
we study the specific example of a false vacuum inflation 
model with a quartic potential || 
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where B is a constant. For suitable B, this model is 
characterized by two separate epochs of inflation. For 
large <p values the false vacuum term is negligible and the 
model becomes a </> 4 potential, while for small values 
the model becomes false vacuum dominated. Depending 
on the value of the parameter B, these two epochs can 
be punctuated by a brief suspension of inflation, while 
the potential becomes temporarily too steep to support 
inflation. Upon reaching the lower part of the potential 
the field fast-rolls until slowed down by friction from the 
expansion. 

The evolution of the quantity 1+e— rj and the comoving 
Hubble wavenumber, aH, is illustrated in Fig. [j] for the 
choice B = 0.55. In this example inflation is suspended 
for about 1 e-fold, but z'/z remains negative for about 
5 e-folds indicating that the field continues to fast-roll 
for some time after inflation restarts. It is during this 
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FIG. 2. The evolution of scalar and tensor perturbations, 
\lZk \ and |Vk|, for a single mode k — 10 -2 which is outside the 
horizon by a factor of 100 when inflation is suspended. The 
solid line is the quantity aH/k, and z'/z < between the two 
vertical dotted lines. The arrow indicates when k = aH, the 
instant of horizon crossing. The absolute normalization of the 
perturbations is arbitrary, though the relative one is correct. 



latter period that scalar modes which have recently left 
the horizon feel the effect of the driving term in Eq. jl^). 

We set the correspondence of scales such that k = 1 
corresponds to the scale which equals the horizon at the 
time when inflation stops. This correspondence is arbi- 
trary, depending on the mechanism for ending inflation 
and in particular the value of <f> at which inflation ends. 

The evolution of perturbations on a particular scale k, 
again for B = 0.55, is shown in Fig. § obtained numer- 
ically using the approach of Ref. E3| . Even though this 
particular mode left the horizon around 7 e-folds before 
z'/z turns negative, the residual lZ' k given by Eq. (|l4|) is 
quickly blown up by the exponentially growing driving 
term. After some time the exponential growth becomes 
important and it drives \7Zk\ through zero and on to a 
much larger amplitude. Although the mode is well out- 
side the horizon at this time, its amplitude is enhanced by 
a factor of around a hundred; we will discuss the physical 
interpretation of this further later. Once z'/z becomes 
positive again normal friction domination resumes, freez- 
ing out \Hk\ at the enhanced amplitude. The standard 
approximations therefore fail by a factor of around a hun- 
dred in this case, which as we will see is by no means the 
worst. The tensor amplitude is unchanged, so the tensor- 
to-scalar ratio is suppressed. 

In Fig. ||we consider two other modes, this time track- 
ing the evolution of |uk|- One mode crosses the horizon 
during the epoch when z'/z < 0. Clearly Eq. (||) does not 
depend explicitly on z'/z, and so the transition from the 
oscillating regime still occurs when k 2 = z"/z. Scalar 
modes that leave the horizon after z'/z turns positive 
again, as is the case in the second mode shown in Fig. |^, 
asymptote to a value that is independent of the influ- 
ence of z'/z turning negative. This is to be expected, 
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FIG. 3. The evolution of two scalar modes, |«k|, as m 
Fig. ^. The arrows again indicate horizon crossing for each 
mode. For the first mode, ki, horizon crossing occurs while 
z' / z < 0, while for the second it occurs after z' / z becomes 
positive again. Both modes asymptote to approximately the 
same amplitude. 



since sub-horizon modes do not feel the influence of any 
background quantities. 

The overall result is that an extremely broad feature 
arises in the scalar power spectrum. This is shown in 
Fig. |i|, where the spectrum for two different choices of 
B has been computed mode by mode, interpolating be- 
tween the two different epochs of inflation. For compari- 
son we plot the amplitude predicted by the slow-roll and 
Stewart-Lyth formulae, where we use the background 
values and slow-roll parameters derived from the numeri- 
cal evolution. For a small range of scales near k — aH | cn( j 
there is some ambiguity as to when the formulae should 
be applied, since these modes leave the horizon twice. We 
plot both possible values. The slow-roll and Stewart- 
Lyth predictions fail by orders of magnitude for many 
e-folds (potentially enough to encompass all scales acces- 
sible to large-scale structure observations in the B = 0.55 
case) , resuming their usual good agreement with the nu- 
merical results soon after z' / z turns positive again. 

An intriguing feature of the spectrum for B = 0.55 
is the very fiat portion on the right, which arises even 
though the disagreement with the slow-roll prediction in- 
dicates that we are nowhere near the slow-roll limit. In 
fact, this is a realization of the more general circumstance 
for a scale-invariant spectrum uncovered by Wands pd| ] 
using duality arguments. During this epoch, the field is 
fast-rolling along a relatively flat section of the potential, 
obeying 
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The solution is 4> oc a , giving z oc r rather than the 
z oc r _1 typical of slow- roll inflation. Nevertheless, this 
relation between z and r leads to a scale-invariant spec- 
trum of scalar perturbations |IT| , just as in the slow-roll 
limit. This decaying solution of 4> is always present but 
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FIG. 4. The scalar power spectrum as determined from 
mode-by-mode integration. The upper panel shows B = 0.55, 
while the lower is for B — 0.3 (note that for inflation to be in- 
terrupted requires B > 0.19 |9|). The scales between the two 
dotted lines correspond to the epoch when z 1 / z < 0. The first 
discontinuity in the analytic expressions occurs in the region 
where modes make multiple horizon crossings, while the sec- 
ond sharp feature in the Stewart-Lyth case occurs due to an 
accidental cancellation of the entire Stewart-Lyth coefficient 
and has no physical significance. 



is usually neglected. Seto et al. Q have shown that in 
extreme cases, when <p becomes equal to zero, the slow- 
roll amplitude of the density perturbation as written in 
Eq. (|l|) will completely fail unless 4> is replaced by the 
corresponding slow-roll velocity <j) s = — V^/SH. 

As the field slows down, one might have expected a 
feature to be produced. However, as long as e is small, 
a necessary and sufficient condition for a scale-invariant 
spectrum is that the entire square-bracketted term in 
Eq. (||) remains constant in time, which taking e <C 1 
reduces to 
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This encompasses both slow- roll (I77I <C 1) and fast- 
roll [r\ > 1) inflation, as well as any smooth transition 
between the two. This condition will be automatically 
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satisfied as long as the inflaton is effectively massless, 
allowing us to neglect the last term of Eq. (Q) . However 
obtaining the flat portion required significant fine-tuning 
of B; notice that with B = 0.3 the potential is much less 
flat in the corresponding region, though still flatter than 
the slow-roll formulae predict. 

In the limit of an instantaneous transition from slow- 
roll to fast-roll behaviour, inflation is no longer suspended 
and we arrive at Starobinsky's model |l3| . In this case the 
inflaton potential is characterized by a sudden gradient 
discontinuity and the power spectrum takes on a similar 
step-like form, but with superimposed oscillations on the 
upper plateau. 

We now return to the physical interpretation of the 
change in the curvature perturbation on super-horizon 
scales. Although ordinarily a single scalar field is as- 
sociated with purely adiabatic perturbations, it can in 
fact support entropy perturbations if its velocity pertur- 
bation does not obey a generalized adiabatic condition 
with respect to the field perturbation (m). Under quite 
general circumstances, however, single field inflation does 
only generate adiabatic perturbations [ p"5| |, with the en- 
tropy perturbation S associated with a general scalar 
field perturbation being non-zero but becoming highly 
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Even in the slow-roll case, care is required in deriving 
the late-time evolution of the entropy, because its sup- 
pression is strong enough that the effect of the 1Z term 
does not become negligible despite its k 2 /a 2 H 2 pref ac- 
tor. If we were able to neglect that term, and taking the 
slow-roll limit, one would find S' /aH ~ —35 implying 
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but we sec that this does fall off faster than 



the last term. The self-consistent solution therefore has 
late-time behavior S ~ e~ 2N , and indeed is what we see 
in our numerical simulations. 

During fast-roll the situation is very different, because 
S is able to grow and so the influence of the final term 
in Eq. ( |25| ) becomes negligible. Making the false vacuum 
(e <C 1) and massless (V^ <C H 2 ) assumptions, Eq. ( p5| ) 
can then be written as 
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suppressed, S 
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, once the mode becomes frozen During fast-roll 77 ~ 3, and the entropy can swiftly grow 



in upon leaving the horizon [14f| . The constancy of the 
curvature perturbation on super-horizon scales in the ab- 
sence of entropy perturbations holds under extremely 
general circumstances 

Our situation is an exception to this. An entropy 
perturbation is the only source of super-Hubble growth 
|^5|,0, and during the phase where z decreases the en- 
tropy perturbation grows until it becomes significant 
enough to source the curvature perturbation. The fur- 
ther above the Hubble length a scale is, the more the 
entropy has been suppressed and hence more time is 
needed for the entropy to become significant, so the ef- 
fect does not extend up to arbitrarily large scales. After 
slow-roll resumes the entropy dies away and the pertur- 
bations become purely adiabatic, remaining so thereafter 
but retaining the shift in amplitude. The influence of the 
scalar field entropy can be studied using the definitions 
of Ref . [03 , where the entropy is given by 
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where A is the perturbation to the metric lapse function 
which is related to the curvature perturbation via the 
constraint equation A — Air^lZ/mp^H 2 0. The en- 
tropy measures the failure of the perturbed velocity to 
match the generalized adiabatic condition with respect 
to the field perturbation itself. It allows Eq. (12) to be 
rewritten as two first-order equations 
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e until it becomes large enough to have a signifi- 
cant impact on 1Z. Once fast-roll ends the k 2 TZ term in 
Eq. (|25| ) is initially small giving the S ~ e~ 3N behaviour, 
but soon after the curvature terms re-asserts itself restor- 
ing the S ~ e~ 2N late-time behaviour. However this 
transition is not of any particular physical significance 
since the curvature perturbation itself has long since ap- 
proached its asymptotic value. 

One might wonder whether the entropy perturbations 
in the scalar field could survive as such after inflation, 
but that does not seem to be possible. Either fast-roll is 
followed by slow-roll as just discussed, or it is followed 
by inflation ending through decay of the inflaton. That 
decay will lead to a purely radiation-dominated universe, 
which is unable to support non-adiabatic perturbations. 
The process of decay of the inflaton would therefore re- 
move the entropy source term. Note also that consider- 
able fine-tuning is required to stay in the fast-roll regime 
for a prolonged period; in order to have inflation the ini- 
tial kinetic energy can at most be of order of the potential 
energy, but then the kinetic energy falls off very quickly, 
and yet for fast-roll to be sustained its contribution to 
the scalar wave equation must remain larger than that 
from the slope of the potential. 



IV. PERTURBATIONS AT THE END OF 
INFLATION 

A second scenario where the standard equations cannot 
be directly applied is at the true end of inflation. This is 
an important regime because such short-scale perturba- 
tions can lead to the formation of primordial black holes, 
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FIG. 5. The scalar power spectrum plotted at various 
stages of its evolution. The numerical spectra are evaluated 
(from top to bottom) at horizon exit, at the end of inflation, 
and at horizon re-entry for each mode. The right-hand edge of 
the plot corresponds to the wavenumber equalling the Hubble 
radius at the end of inflation. The Stewart-Lyth and slow-roll 
spectra are evaluated, as usual, at horizon exit. 

which gives the most important constraint on the late 
stages of inflation. There are two standard mechanisms 
for ending inflation, one being the hybrid method of an 
instability in another direction in field space and the sec- 
ond being the breakdown of slow-roll (see Ref. || for a 
review) . A particularly relevant case is in hybrid inflation 
models such as the running-mass model [l6| ], where the 
slow-roll approximation is only well respected over a lim- 
ited range of scales. In the running-mass model slow-roll 
inflation comes to an end due to -q growing, but inflation 
may continue in the fast-roll regime until an instability 
is reached. 

Typically the standard formulae for the perturbations 
will break down near the end of inflation, because of a 



failure to reach the asymptotic limit (see e.g. Ref. (17|) 
and often because the slow-roll approximation is not ac- 
curate. For example, in the running-mass model Eq. dl"" 
becomes 



(27) 



and so r\ = 1 is the last point at which constraints based 
on the shape and size of a blue inflationary power spec- 
trum can be reliably applied without resorting to mode- 
by-mode integration of the power spectrum |1| . 

We are not able to study the hybrid case due to 
the complexity of the multi-field dynamics (though see 
Ref. [jl4| for a general formalism for doing so), so we re- 
strict our attention to single-field models ending by vio- 
lation of slow-roll. This is not in fact the most interesting 
case as typically such models have a red spectrum where 
small-scale perturbations are not very significant, but il- 
lustrates the main physics. For simplicity we study the 
quadratic chaotic inflation model V(<f>) oc 4> 2 . 

In Fig. |^ we plot the scalar power spectrum for the last 
few e-folds of modes to leave the horizon before inflation 



ends. Because the asymptotic regime is not adequately 
reached at any stage, there is no preferred choice as to 
when to plot the amplitude, and we evaluate it at three 
stages: horizon crossing during inflation, at the end of 
inflation, and at horizon re-entry after inflation once the 
background field begins oscillating. Notice that the am- 
plitude at horizon exit is typically much greater than the 
value at the end of inflation, and that the slow-roll and 
Stewart-Lyth approximations indeed resemble the latter 
rather than the former. 

We see that the slow-roll and Stewart-Lyth formulae 
underpredict the asymptotic and re-entry amplitudes of 
perturbations which exited the horizon close to the end 
of inflation. One therefore expects an enhancement of 
primordial black hole production, though a detailed cal- 
culation would need to track the derivative of the cur- 
vature perturbation at re-entry as well as its amplitude. 
This result suggests that normally quoted constraints are 
on the conservative side, though typically the correction 
would not be large. We mention additionally that for 
perturbations which do not reach the asymptotic regime 
there are questions as to how the quantum-to-classical 
transition might take place (see e.g. Ref. jl9|]); we will 
not however attempt to address this here. 



V. DISCUSSION 

The accuracy of the usual analytic expressions for the 
inflationary power spectrum depends on scales evolving 
smoothly through horizon crossing and into the asymp- 
totic regime k 2 <C z"/z. In this paper, we have investi- 
gated two situations where this is not achieved, one being 
a temporary end to inflation and the other the true end. 
In the former case we have seen that the modes can have 
a very complicated evolution, including the possibility 
of amplification on super-horizon scales via an exponen- 
tial driving term corresponding to an entropy perturba- 
tion in the scalar field. Such behaviour can be traced 
to the complicated evolution of the scalar pump field z, 
which no longer grows monotonically. The net effect is 
the insertion of a broad feature into the power spectrum 
that can only be computed by mode-by-mode integra- 
tion, and which can differ wildly from the slow-roll and 
Stewart-Lyth approximations. Typically the features are 
sufficiently non-scale-invariant to be excluded already on 
astrophysical scales, but we have also seen that a very 
flat spectrum can be obtained while far from slow-roll, 
confirming an analytic analysis by Wands JlT[ . 

The second scenario we studied was the true end of in- 
flation, where the failure to reach an asymptotic regime 
means that perturbations re-enter the horizon with a 
higher amplitude. The most interesting physical con- 
sequence of such modes is in primordial black hole for- 
mation, and this result indicates that earlier treatments 
assuming the slow-roll formulae are somewhat conserva- 
tive, as the mode equation solutions indicate that these 
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formulae underestimate the perturbation amplitude at 



horizon re-entry. 
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